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ON  THE  PROPAGATION  OF  THE  PLASTIC  DEFORIiATION 


PRODUCED  BY  AN  EXPA!ffiING  CYLINDER 


Abstract 

In  the  present  paper  it  is  assumed,  in  accord  mth  Bethe's  model 
for  armor  penetration,  that  the  plastic  deformation  produced  during 
armor  penetration  is  similar  to  that  produced  by  an  expanding  cylinder. 
The  paper  deals  with  the  deformation  produced  in  a  plate  which  is  thick 
enough  so  that  the  plate  does  not  get  appreciably  thicker  near  the  ex¬ 
panding  cylindrical  hole.  In  the  first  portion  of  the  paper  the  rela¬ 
tions  between  the  stresses  and  the  strains  are  considered  and  the  wave 
equations  which  govern  the  motion  of  the  material  virhen  it  is  rapidly 
deformed  are  derived.  It  is  found  that  in  the  case  of  the  thick  plate 
an  elastic  wave  diverges  radially  from  the  expanding  cylinder,  that 
this  is  followed  by  a  plastic  v/ave  and  that  the  elastic  and  plastic 
wave  velocities  do  not  differ  very  much.  In  the  second  portion  of  the 
paper  an  approximate  expression  for  the  displacement  is  obtained  for  the 
case  where  the  deformation  is  elastic.  In  the  last  section  of  the  paper 
a  method  of  numerically  integrating  the  wave  equations  is  given.  The 
method  is  applied  to  a  particular  numerical  example,  and  the  displace¬ 
ment  and  the  shearing  strain  produced  by  a  uniformly  expanding  cylinder 
are  calculated  and  plotted  at  various  times  during  the  expansion.  The 
calculation  shows  that  the  shearing  strain  in  both  the  elastic  and  the 
plastic  regions  increases  as  we  go  towards  smaller  radii.  A  discon¬ 
tinuity  in  the  shearing  strain  is  found  at  the  boundary  between  the 
elastic  and  plastic  regions.  The  calculation  indicates  that  the  com¬ 
pressibility  and  the  density  play  an  important  role  in  determining  the 
magnitude  of  the  stresses  developed  in  a  thick  ..late.  A  program  for 
future  investigation  is  suggested. 


1 .  Introduction 
1  / 

Bethe—  has  calculated  the  stresses  and  strains  produced  in  armor 
plate  by  a  projectile  under  the  assumption  that  the  actual  three-dimen¬ 
sional  problem  can  be  approximated  by  a  treatment  involving  only  two 
dimensions.  Diu-ing  penetration,  a  pointed  shell  will  produce  large  dis¬ 
placements  of  the  surrounding  material  away  from  the  axis  of  the  shell. 


y  All  numerical  references  are  to  the  list  of  references  given 
at  the  end  of  this  report. 
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Bethe  assumes  that  in  a  thick  plate  the  displacements  parallel  to 
the  axis  of  the  shell  are  negligible  compared  with  displacements 
perpendicular  to  the  axis.  The  stresses  and  strains  are  calculated' 
by  assuming  that  the  projectile  produces  a  cylindrical  hole  in  the 
plate.  As  the  penetration  of  the  shell  proceeds,  the  hole  mdens. 

The  initial  radius  of  the  hole  is  zero  and  its  final  radius  is  the 
radius  of  the  shell. 

The  forces  that  displace  material  away  from  the  axis  of  the 

shell  must  accelerate  the  material,  and  they  must  overcome  the 

stresses  produced  by  the  deformation.  Bethe  recognizes  that  iper- 

tial  forces  will  play  a  role  in  the  problem,  but  he  neglects  them 

in  his  treatment.  '  Calculations  on  the  propagation  of  plastic  de- 

2/ 

formation  in  steel  wipes—  yield  kinetic  energies  that  in  all  cases 
exceed  hO  percent  of  the  total  energy.  It  is  to  be  expected  that 
the  kinetic  energy  may  represent  a  smaller  fraction  of  the  total 
energy  in  the  case  of  cylindrical  waves,  but  one  would  not  expect 
the  kinetic  energy  to  be  entirely  negligible. 

In  his  calculation  Bethe  has  given  a  static  treatment  of  the 
problem.  This  means  that  he  calculates  the  stresses  and  the  strains 
present  in.  the  plate  at  an  infinite  time  after  the  penetration  of 
the  projectile.  Since  the  initially  intense  disturbance  around 
the  hole  eventually  spreads  oyer  a  larger  area,  it  is  clear  that 
the  maximm  stresses  and  strains  produced  in  the  armor  are  under¬ 
estimated  in  Bethe 's  calculation. 

G.  I.  Taylor  hap  also  considered  the  problem  of  armor  penetra¬ 
tion.^  His  first  paper  gives  a  static  discussion  of  the  deformation 
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produced  in  a  thin  plate  by  an  expanding  cylinder.  Taylor's  relations 
connecting  the  stresses  and  the  strains  differ  from  those  used  by  Bethe, 
In  his  second  paper  Taylor  gives  a  dynamic  treatment  of  the  enlarge¬ 
ment  of  a  hole  in  a  thin  plate  at  high  speeds]  that  is,  in  this  paper 
he  includes  inertial  forces.  The  assumption  is  made  that  the  tensile 
stress  which  acts  in  the  plane  of  the  plate  and  perpendicular  to  the 
radial  direction  is  zero. 

In  the  present  paper  we  shall  discuss  the  elastic  and  plastic 
Yraves  set  up  in  a  thick  plate  by  an  expanding  cylinder.  We  shall  ob¬ 
tain  first  the  necessary  viave  equations.  An  approximate  expression  for 
the  radial  displacement  in  the  elastic  case  will  then  be  found.  Next, 
a  method  for  numerically  integrating  the  mve  equations  will  be  de¬ 
scribed  and,  finally,  some  results  obtained  by  nuraerical  calculation 
vdll  be  given. 


2.  The  wave  equations 

From  the  symmetry  of  the  problem  we  see  that  the  displacement  pro¬ 
duced  by  an  expanding  cylinder  in  a  thick  plate  is  entirely  radial.  The 
equation  of  equilibrium  can  then  be  wrritten,— 


^<ir 


«r  -  ^ 


(1) 


where  u  is  the  displacement,  ^  is  the  density  of  the  material,  r  and  G 
are  plane  polar  coordinates,  and  cTg  are  tensile  stresses  in  the 
directions  r  and  G,  and  t  is  the  time.  If  we  assume  that  the  displace¬ 


ment  is  small  compared  mth  r  then,  according  to  Timoshenko 
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we  have 


-h- 


p 

~  ^r©  ”  ^rz  ~  ^©z  ”  ^ 


(2) 


'Here  Ej,,  Cg  and  6^  are  the  tensile  strains  in  the  directions,  r,  © 
and  z,  respectively!  7j^q)  7-^2.  ^©z  shearing  strains! 

and  the  z-axis  is  the  axis  of  the  shell.  Equations  (2)  hold  only 
for  the  case  where  the  displacement  is  completely  radial.  It  might 
be  well  to  point  out  that  actually  u  is  not  small  compared  mth  r 
near  the  axis  of  the  shell  —  that  is,  near  the  z  axis.  Mo  treat¬ 
ment  has  as  yet  been  given  that  deals  correctly  mth  this  region. 
We  shall  use  Eqs.  (2). 

\ 

(a;  Elastic  deformation.  —  A  mve  equation  in  the  variable  u 
can  be  obtained  if  we  assume  relations  between  the  stresses  and 
the  strains.  If  the  deformation  is  elastic,  the  relations  are 
given  by  Hooke's  law.  Thus, 


(3) 

*  (T-  =  2G  ,  -  (T  +  +  (T  =  3K(6  +  £, ) , 

re  r'^  r©z  r© 

where  G  is  the  modulus  of  rigidity  and  K  is  the  modulus  of  volume 
expansion.  Note  that  does  not  appear  in  Eqs.  (3)!  in  the  parti¬ 
cular  problem  considered  here,  is  zero.  If  vre  use  Eqs.  (2)  and 
substitute  values  obtained  for  the  stresses  from  Eqs.  (3)  into 
Eq.  (1),  we  obtain  the  mve  equation, 

3^u  /3K  +  Ug  \ /3^u  1  3u  u\ 

3/0 
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The  velocity  of  propagation  of  these  elastic  v/aves  is 

Cl  =  V(3K  +  kG)/3p.  (5) 

If  we  calculate  c^  for  steel  for  the  case  where  G  =3.140x10^^  Ib/ft  sec^, 
K  =  I.OI46  X  10^^  Ib/ft  sec^  and  .I18O.7  Ib/ft^,  we  find  that 
Cl  =  1  9^1  70  ft/sec. 

6/ 

According  to  Mohr's  criterion—  the  material  at  a  specified  point 
mil  experience  plastic  deformation  if  the  largest  component  of  the 
shearing  stress  at  the  point  exceeds  the  limit  g-Aj,  where  A  is  the 
yield  stress  obtained  in  tensile  tests  made  on  long  thin  mres.  Huber, 

7/ 

Henckey  and  von  Mises-  postulate  that  plastic  deformation  of  a  poly¬ 
crystalline  material  begins  when  the  principal  stresses  become  large 
enough  to  satisfy  the  relation. 


(T,  -  (T^)^  +  ((^3  -  ^3)^  +  (<r3  -  CTi)^  =  2\"  .  (6) 


The  conditions  of  Mohr  and  von  Mises  do  not  differ  very  much.  The  Tforst 
agreement  is  obtained  in  the  case  of  pure  shear.  In  this  case  the 
largest  component  of  the  shearing  stress  found  by  using  Mohr's  criterion 
is  13.3  percent  lower  than  the  value  obtained  by  using  the  von  Mises 
condition.  Yfe  shall  use  Mohr's  condition  to  determine  whether  or  not 
the  material  at  a  given  point  is  experiencing  plastic  deformation. 

(b)  Plastic  deformation.  —  We  shall  assume  that  the  follomng 
relations  are  valid  for  material  that  is  being  plastically  deformed; 


-  _  1 

,  _  ) 

■q"  D  L”g 


rp]  a; 


r  -  a((r  +  (T 


0  “o, 

4.  (f 


f  3; 

2  D 


.  (1. 

rV  - 

5'!| 

(P  4-  op. 

\E 

D/  Zq 

Is 

0/1 

L  ^’0  QoJ 

(7) 
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In  Eqs.  (7),  QTy,  }  (Tq  and  r  are  the  stresses  required  to  initiate 
.plastic  deformation.  The  constants  E  and  v  are  Young's  modulus  and 
the  Poisson  ratio.  The  constants  D  and  ^  describe  the  plastic  be¬ 
havior  of  the  materialj  the  7ra.y  in  which  they  are  determined  from 
experimental  data  1/0.11  become  clear  as  ive  proceed.  ¥e  shall  wait 
until  later  to  discuss  the  validity  and  meaning  of  these  equations. 
Adding  Eqs.  (7)  we  obtain 


1-  2o< 


^r  ^  +  ‘-^z - D — ^  ‘^zi 


1  -  2v  1  -  2a'^ 
E  “  D  i 


[(Tj.  +cr  ].  (8) 

”o  ^0 


In  the  elastic  cage. 


^  ^  +  ^z  '  3T  C<^r  E^r  ^  “"z^^  (9) 

where  1  /K  is  the  compr8S.sibility.  life  shall  assume  that  Eq.  (9)  is 

8  / 

also  valid  for  a  material  which  is  being  plastically  deformed. - 
Since  Eq.  (8)  must  be  the  same  as  Eq.  (9),  we  have 

This  equation  can  be  used  to  determine  a.  once  D  is  known.  D  can  be 
obtained  from  the  stress-strain  curve  of  the  material  in  tension. 
Equations  (7)  assume  that  thfs  sl^ress  strain  curve  is  of  the  form 
shown  in  Fig.  1.  The  slope  of  the  initial  elastic  portion  of  the 
curve  is  E.  The  slope  of  the  stress  strain-curve  in  the  plastic 
region  above  the  knee  is  Equations  (7)  are,  of  course,  valid 

only  if  r^.,  Iq  and  (Tg  principal  stresses. 


r 


Fig.  1.  The  form  of  the  stress-strain  curve 
assurued.  in  Eqs.  (7). 


Ros  and  Eichinger,  Nadai,  and  others—  have  shoYm  that  relations 
of  the  type  of  Eqs.  (7)  are  in  accord  ivith  experimental  tests  made  on 
polycrystalline  materials.  The  relations  are  so  constructed  that  the 
directions  of  the  principal  strains  coincide  vlth  the  directions  of 
the  principal  stresses.  In  addition,  the  equations  are  so  arranged 
that  the  figure  consisting  of  Mohr’s  three  principal  strain  circles 
remains  continuously  similar  geometrically  to  the  figure  made  up  of 

9/ 

the  three  principal  stress  circles  (see  Nadai).~  Our  relations  are 
someidiat  more  general  than  theirs  since  they  took  K  to  be  infinite  — 
Otis  then  0.5  by  Eq,  (10).  The  fact  that  K  is  not  infinite  is  of 
importance  in  the  problems  treated  in  the  present  report. 

Let  us  consider  the  plastic  deformation  vihich  occurs  in  a  thick 
plate.  For  this  case  we  have: 


=  v((rr^  +  re^),  (n) 

Equations  (3)  hold  at  the  boundary  between  the  elastic  and  plastic 
regions  and,  since  we  still  assume  thab  the  displacement  is  radial  in 
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the  plastic  region^ Eq,  (2)  is  still  valid  there.  Using  the  last 
of  Eqs.  (2)  we  can  write  Eqs.  (?)  as  .  * 

«'e  -  «-r  -sftEe  -  ti-l  +  2(g  -F)[E.e^  - 

-  iTg  --2F&Q  -  2(G  -  F)Eg  ,  '  (12) 

G-j.  -  (Tg,  =  +2F£^  +  2(G  - 

where  En  and  £„  are  the  strains  present  when  the  material  begins 
^  0 

to  deform  plastically,  and  where 


and 


F  = 


D 

2(1  +  oc) 


°  “  2(1  +  V)> 


G  being  the  elastic  shear  modulus.  In  addition  we  have  from.Eq.  (9), 
since  =  0, 


=  3K[£e  +  6^^]. 


(13) 


The  wave  equation  for  plastic  waves  is  obtained  by  using  Eqs.  (1), 
(2),  (12)  and  (13)  in  the  same  way  that  Eqs.  (l  ),  (2)  and  (3)  were 
used  to  obtain  the  elastic  wave  equation.  The  wave  equation  which 
results  is 


a^u  nK  +kF\/'fM  1  '3u  u\  2(G  -  F)y^ 

c>t^  \  3f-^  /  n  'ftn  I*  /  /=>r 

where  7-  ~  Sr,  “  '  1^1^®  plastic  mve  velocity  is 

w  ^  O 


(1U) 


Cg  -  V(3K  'y  hF)/3fi 


(15) 
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Using  F/G  =  ^/6.9h  foi”  steel,  ive  find  that  Cg  =  ^S,h^O  ft/ssc.  The 
ratio  C;j/c]_  is  therefore  0.80.  For  steel  of  the  quality  used  in 

armor  plate  A  =  I4.6I  x  10®  Ib/ft  sec^,  =  1|.27  x  10“^  and 
2>o(G  -  F)//o  =  8.20  X  105  ft^/sec^ 

(c)  Plastic  deformation  when  the  slope  of  the  stress-strain 
curve  is  some  function  of  the  largest  shearing  strain.  —  Yfe  can 
generalize  our  vra.ve  equation  by  assuming  that  the  material  folloTfs 
a  stress-strain  curve  which  can  be  introduced  into  the  equations  as 
an  arbitrary  function.  If  y/e  suppose  that  the  slope  of  our  stress 
strain  curve  is  s  one  function  of  the  largest  shearing  strain,  then 
the  eqmtions  relating  the  stresses  and  the  strains  may  be  y/ritten 
in  a  form  analogous  to  Sqs.  (3)  as 

^  '^^^9  ~  ^r^  ”  ^r^^ 

(16) 

S  *  \  ‘•el- 

Ros  and  Eichinger  have  sho™  experimentally  that  the  D  and  F  of 
our  Eqs.  (7)and  (12)  actually  are  not  constant  but  depend  somewhat 
on  the  variable  quantity,  [(Eq  -  -  6^)^  -f  (£^  - 

Yfe  shall  assume  that  only  the  largest  shearing  stress  plays  an  im¬ 
portant  role.  This  assumption  is  related  to  the  results  of  Ros  and 
Eichinger  in  the  same  manner  as  Mohr's  condition  for  the  onset  of 
plastic  deformation  is  to  the  Huber,  Honckey  and  von  Mises  condition. 
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The  function  f  may  be  determined  from  static  experi¬ 

ments  in  the  following  way."^  In  general,  the  equations  relating 
the  principal  strains  and  the  principal  stresses  can  be  written, 

n'-  -  «3), 

-  ^3  -  2(5^  -  -  £3), 

“  2(6-3  —  ^i)f(^l  “  ^3) 

and 

'^1  ^2  “  3K(6-j_  +  £3  ^3)^ 

where  ~  ^3  is  the  largest  shearing  strain.  These  relations  will 
then  be  valid  in  the  case  of  a  tensile  test  on  a  long  thin  ware  ex¬ 
tending  in  the  Xi-direction.  In  this  case, 

2  =  ^3=0  .  ' 

and,  hence,  by  the  second  of  our  general  equations, 

~  ^3* 

The  equations  relating  the  principal  strains  and  principal  stresses 
become 

^  2(e^  -  e^)f(e^  -£3),  =  3K(&i  +  2  £-3). 

Eliminating  £3  from  these  equations,  wfe  get 


This  method  was  suggested  by  M.  P.  IS/hite  and  H,  F.  Bohnenblust. 
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Thus  by  measuring  ^  tensile  test  and  by  plotting  the 

function  on  the  right  as  a  function  of  {^^x~  function  f 

can  be  found,  If  K  can  be  taken  as  infinite  —  that  is,  if  the  sub¬ 
stance  is  incompres sable  —  then  the  foregoing  equation  becanes. 


Equations  (16)  are  so  constructed  that  the  directions  of  the  principal 
strains  coincide  with  the  directions,  of  the  principal  stresses.  In 
addition,  the  equations  are  so  arranged  that  the  figure  consisting  of 
Mohr's  three  principal  strain  circles  remains  continuously  similar 

geometrically  to  the  figure  made  up  of  the  three  principal  stress 

.  .  9/ 

circles  (see  Nadai}.~ 

Proceeding  as  before,  using  Eqs,  (1),  (2)  and  (16),  we  find  that 
the  wave  equation  for  a  material  in  which  the  slope  of  the  stress- 
strain  curve  depends  on  the  largest  shearing  strain  is 


^^u  [K+^f-|(2Cr-2e)f’ I  K  +  ^f+|(2£j.-  ]|~1  ^u 


-f  +  -( 

3  3 


-—--2  .  (17) 
r  ar  r 


In  this  equation  f '  represents  the  deri-vative  of  the  function  f  mth 
respect  to  its  argument  Equation  (17)  is  complicated  be¬ 

cause  the  coefficients  of  the  partial  derivatives  depend  on  u.  Von 

io/ 

Karman — '  has  discussed  the  propagation  gf  plastic  deformation  in 


wires  when  the  stress-strain  curve  is  given  as  a  f'onction  to  be  deter¬ 
mined  by  experiment,  'fe  have  nqt  been  able  to  apply  von  Karman 's 
method  of  integration  to  Eq,  (I?)* 
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(^)  Plastic  deformation  in  a  region  Tiyfaere  the  shearing  strain 
decreases  -with  the  passage  of  time.  —  After  the  conical  head  of  the 
projectile  has  passed  ccmpletely  through  the  armor  plate,  the'  stresses 
in  the  ■vicinity  of  the  cylindrical  hole  will  begin  to  decrease.  We 
shall  next  obtain  a  vrave  equation  appropriate  for  this  region.  If  a 
metal  is  plastically  defomed,  moving  from  A  to  B  to  C  on  the  stress- 
strain  curve  shown  in  Fig,  2,  and  if  the  largest  shearing  stress  is 


Fig.  2.  The  stress-strain  curve  for  plastically  de¬ 
formed  material  including  a  region  ®  when  the  shearing 
stress  is  decreasing.  Note  that  CD~Ts  parallel  to  the 
elastic  portion  of  the  curve 

then  decreased,  the  shearing  strain  also  decreases 5  the  stress  and 
strain  move  down  along  the  straight  line  OT.  The  line  CD  is  paral¬ 
lel  to  the  elastic  portion  of  the  stress-strain  curve |  that  is,  CD 
is  parallel  to  If  and  are  the  strains  at  point  C  of 

the  curve  and  if'  7x  is  the  shearing  strain  at  that  point,  then  the 
equations  relating  stress  and  strain  along  the  line  CD  are 

(rQ-(r^=2G(£e-£^)-2(G-F)(7i-?'o),  e-^-6-Q=«2G£e+2(G-F)(Se^-£Q^), 

(r^-(r2=  +2G£^-2(G-F)  3K  £9). 


(18) 
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The  -wave  equation  obtained  by  using  these  equations  is 

9^u  _  /3Kj^y^  ^  2(G  -  F)a.i  -7o) 

3t^  \  3p  r  or  r  /  ^r 

The  mve  velocity  is  just  the  velocity  of  the  elastic  mves. 

3.  The  elastic  case 

In  order  to  obtain  seme  idea  of  the  nature  of  the  functions  that 
appear  in  the  solutions  of  the  v;ave  equE.tions,  we  have  attempted  to 
solve  a  problem  involving  only  elastic  vra-ves.  Consider  a  thick  plate 
of  infinite  extent,  through  the  center  of  v^hich  passes  a  cylindrical 
hole  of  radius  b.  We  shall  assume  that  external  forces  act  against 
the  surface  of  the  hole  so  that  its  radius  increases  with  the  passage 
of  time.  The  radial  displacement  at  b  is  supposed  to  vary  vdth  time 
in  the  following  fashion: 

u(b)  =  0,  for  t  <  Oj 

u(b)  =  Vgt,  for  0  <  t  <  T|  (20) 

u(b)  =  VqT,  for  T  <  t. 


where  v^  and  T  are  constants.  Since  the  deformations  are  taken  to  be 
elastic,  the  w^ave  equation  is 


^^u  _  ^2  ^  1_  ^  u  ~j 

3t^  _3r^  r  3r 


(U) 


Bathe —  has  given  a  theoretical  discussion  qf  a  very  similar  problem. 


¥e  have  extended  his  treatment  to  higher  ordeps  of  approximation,  and 


we  have  used  the  method  to  obtain  an  expressj.on  for  u  that  is  valid 
in  the  portion  of  the  vrave  sent  out  while  the  cylinder  is  expanding. 


life  find  that 


n  Vob 


b  fr  / 

- fl+— - 


^j,6  \  b  /  b^  b^  $b^ 

- + - )-*•  -1-  —  1 1  +- - h - ^ — + - 

I.  1.  (i.  ^  /  \  o  o  1,  I.  A 


- / 1  + - (- - + - )-'•  ~  j  1  +“ - ^ ^ — + - -  +  ' 

2cit  b  \  Iic^t^  Sc^t^  6Ucit*^  ''  r  V  hc^t^  8cit^  6Ucit'^ 


b^r  log  g 
2G^t^ 


3(r^  +  b^)  l5(r^  +  3r^b^+ b^) 


TTb^r  3(r^  -  b^)  ^  5(r^  -  b^)  ^ 
c^t^  .  h Cit^  ^  2  c^t^ 


/If  f  \f  2CT.t 

+2/ - )(  Z'(3)  -log  — 

\rcj^t  be  it/  \  b  , 


c^t  yre^t^  bc^t^  /\  b 


15/  5b^  (br^+  3b^r^  +  b^r) 


k  Vre^t^ 


boft^ 


•)lZ'(7)-log  - 

b  / 


In  this  equation  Z‘(a)  =  vdiere  I  (a)  is  the  gamma  func- 

da 

tion  evaluated  for  argument  a.  The  equation  is  valid  if 
Gi(t  -  T)  <  r  <  Cit.  The  eemplexity  of  Eq.  ( 21 )  indioates  that  it 
would  be  diffiGult  to  find  an -analytical  expression  for  the  dis¬ 
placement  in  the  case  where  both  elastic  and  plastic  waves  are 


present.  The  sitijation  would  then  be  more  complicated  than  it  is 


-In¬ 


here  since  the  solutions  obtained  in  the  elastic  and  plastic  regions 
must  be  fitted  together  at  the  boundary  between  the  two  regions. 


k.  Numerical  calculations 

In  our  opinion  a  qualitative  understanding  of  the  waves  produced 
by  an  expanding  cylinder  is  most  easily  obtained  by  a  numerical  inte¬ 
gration  of  the  wave  equations.  In  this  section  a  method  for  numeri¬ 
cally  integrating  the  wave  equations  is  described,  and  some  results 
obtained  by  numerical  oalciijLation  are  presented. 

Let  us  assume  that  the  cylindrical  hole  is  expanded  rapidly  enough 
to  cause  plastic  deformatipn.  In  this  case  both  an  elastic  and  a 
plastic  Tjave  are  produced.  The  elastic  wave  runs  out  ahead  .of  the 
plastic  wave  Just  as  it  does  when  a  large  force  tending  to  lengthen 
a  long  wire  is  suddenly  applied  Ibo  one  end  of  the  ifnre.  The  wave 
equations  wrhich  must  be  considered  are: 


5t' 


2 

Cl 


it^ 


Co 


'fn 


1  Ou 
r  c'r 


+  i  ^ 

r  3r 


u 


Ul  __  -  F) 


2  I 

r  J 


(Elastic)  (h) 
(Plastic)  (lh) 


These  tvro  differential  equations  can  be  written  approximately  as  finite 
difference  equations.  To  do  this  we  shall  consider  only  the  displace¬ 
ments  at  regularly  spaced  instants  of  time  and  at  regularly  spaced 
radii.  The  time  interval  will  be  called  A|  the  space  interval  in  the 
elastic  case  we  shall  take  to  be  c^Aj  and  the  space  interval  in  the 
plastic  case  mil  be  C2A.  ?|fe  shall  denote  by  Uj^j  the  displacement 
at  the  i  th  instant  of  time  and  at  the  radius  r^.  In  this  notation 
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the  differential  equation  for  the  elastic  case  becomes 


^(i+1)  j“2uj_  j+u(  j__i ) 

^i(j+1)  j-1  )^^  j 

“Ij 

A^  ^ 

L  c|a2  r. 

V  '2  c^A  /  r^ 

J 

Solving  this  equation  for  ^  we  find  for  the  elastic  case, 

2 

c^A/-  »  (ciA)  Uj_^ 

^(i+1)j''"^(i-1)j’*'^i(j+lf^i(j-l)'*'  )“^i(  j-1  )y  •  (22) 


Similarly  in  the  plastic  region  vre  find 

c,A/ 

^(i+1 )  f -^(i-1 )  )  "^i  ( j-1 ) ^ 

(23) 

where  P  =  2}'q(G  -  F)/*^.  The  method  of  calculation  can  be  under¬ 
stood  mth  the  aid  of  Fig,  3*  In  this  figure  the  region  that  lies 
above  the  solid  line  represents  the  elastic  region;  the  region  below 
this  line  represents  the  plastic  region. 

Since  the  maximum  velocity  of  any  wave  in  the  material  is  C]_ 
and  since  the  material  is  undist\irbed  until  the  time  t  =  0,  it  is 
evident  that  the  displacements  in  the  region  above  the  line 
r  “  b  =  c^t  are  zero.  We  shall  assume  that 


(cgA)  PA 

V  A 


u(b)  =  0,  for  t  <  0; 
u(bj  =  Vgt,  for  t  >  0. 


(2U) 


Thus,  the  displacements  -Uj^q  are  also  knovi'-n.  According  to  Eqs.  (22) 
and  (23)  a  particular  displacement  can  be  calculated  if  four  neigh¬ 
boring  displacements  are  known.  The  position  of  these  displacements 
on  the  diagram  is  simply  related  to  the  position  of  the  displacement 


'T/me 
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being  calculated.  For  example,  if  one  wishes  to  calculate 
(see  Fig.  3),  the  displacements  required  lie  in  an  inverted  T 
just  above  shearing  strain  at  any  point  i,^  is  given 

by  .  ' 

.  -ii  _ -iiJj-U - iitli  .  (25) 

There  are  several  ways  in  which  the  calcjalation  may  be  per¬ 
formed  j  we  have  found, the  displacements  as  follows.  The  elastic 
displacements  are  calculated  first.  The  value  of  u^  is  found  by 
requiring  tl^at  the  average  shearing  strain  between  the  point  y 
and  the  pbint  A  of  our  diagram  must  be  equal  to  the  shearing  strain 
which  is  required  to  produce  plastic  deformation.  Thus> 

Uii  Uii  -  0 

Tb+Acs  A(ci  -  Cs)  * 

The  displacement  u^g  then  be  calculated  at  once  because,  of 
the  four  displacements  necessary  for  its  calculation,  only,  Un 
differs  from  zero.  Knowing  Ugg  we  can  calculate  U33,  and  so  on. 
Since  the  elastic  wave  travels  faster  than  the  plastic  wave,  it 
soon  becomes  necessary  to  introduce  another  diagonal  row  of  dis¬ 
placements  in  the  elastic  region.  We  have  adopted  the  practice 
of  inserting  an  additiona.1  row  whenever  the  difference  in  radius 
between  the  smallest  elastic  radius  at  a  particular  time  differs 
frcm  the  largest  plastic  radius  by  more  than  A(cj  +  Cg).  Sup¬ 
pose  that  the  displacement  at  the  top  of  a  new  diagonal  row 
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is  Then  is  found  by  making  the  average  shearing  strain  be¬ 

tween  the  points  (a,b)  and  (a,b+  1)  equal  to For  instance,  u^^ 
in  the  elastic  region  is  given  by 

,  ^65  ~  ^66  ^65  ^  ^66 

V  o  .  • 

Aci  r^  +  r^ 


The  other  members  of  the  nev^  diagonal  rov:  and  so  forth,  are 

then  calculated  successively  using  the  finite  difference  equation, 

Eq.  (22) . 

After  the  displacements  in  the  elastic  region  have  been  evalua¬ 
ted  for  times  between  t  =  0  and  any  given  later  time,  the  displace¬ 
ments  in  the  plastic  region  can  be  calculated  for  the  same  range  of 
time.  Consider  first  the  displacements  along  the  line  r-b  =  C2(t-A). 
In  this  diagonal  roi,?,  U21  can  be  found  immediately  in  terms  of 

^io(plastic)  ^ii(elastic)*  displacement  cannot  be  evalu¬ 

ated  at  once  because  the  elastic  displacement  u^  is  not  knoi^n  (see 
Fig.  3)-  The  same  situation  occurs  all  along  this  diagonal  row.  For 
example,  suppose  that  wo  have  found  u^^^  and  now  aush  to  evaluate 

the  displacement  Uo„  can  bo  found  only  if  wo  know  Ur,,  and  u,^. 

c>  7  D  Ji  r 

ITo  have  evalu.tod  the  displacements  at  points  such  as  B,  D,  E  and  F 
by  fitting  the  nearest  elastic  displacements  to  a  linear  relation 
between  displacement  and  radius.  The  linear  relation  found  in  this 
way  is  then  used  to  calculate  the  desired  displacements.  For  example, 
it  is  assmed  that  the  displacement  increases  linearly  mth  distance 
as  we  go  from  point  C  to  point  2,2,  The  displacements  u^  =  0  and  U22 
are  used  to  determine  the  constants  in  the  linear  relation  and  u^  is 

D 
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then  calculated  using  the  linear  relation.  Similarly  u^y  and  Uy^ 
determine  the  linear  equation  that  is  used  to  find  Up  and  Ug. 

The  displacements  in  additional  diagonal  rows  lying  successively 
farther  frcm  the  first  diagonal  plastic  row  are  then  calculated 
\xntll  all  of  the  plastic  displacements  bet’.veen  t  =  0  and  the  pre¬ 
viously  specified  time  are  knov/n. 

Since  the  energy  density  associated  ivith  any  particular  place 
in  our  TC.ves  decreases  as  they  spread  outY.Ci.rds,  it  is  not  sur¬ 
prising  to  find  that  the  shearing  strain  at  the  outer  edge  of  the 
plastic  wave  decreases  Yvdth  time.  Thus,  if  we  evaluate  the  average 
shearing  strain  over  the  interval  between  the  smallest  elastic 
radius  and  the  largest  plastic  radius,  this  strain  is  for  small 
times  larger  than  but  eventually  decreases  until  it  equals 
l^Tien  this  shearing  strain  bocemes  smaller  than  the  first 
diagonal  plastic  rovr  must  be  replaced  by  a.  diagonal  elastic  roYir. 

ITe  shall  show  hovY  the  replacement  is  made  using  the  diagram.  Let 
us  suppose  that  y«3  have  calculated  and  . 

In  our  calculs-tions  in  the  elastic  region  wo  ha.ve  started  a  new 
diagonal  elastic  roT,v  at  ^  since  there  the  difference 

betYveen  the  smallest  elastic  radius  r.jQ  ^he  largest 

plastic  radius  r.jQ  9(plastic)  become  greater  than  A(ci  +  Cg). 

We  shall  suppose  that  further  calculation  shOYvs  that  the  shearing 
strain  between  (8, 7) (^i^^gtic)  (plastic)  Sweater  than 

3'o,  Yfhile  the  shearing  strain  betYveen  and  8) (plastic) 

is  less  than  We  therefore  discard  our  value  for  '^9S(piastic) 

and  move  the  edge  of  the  plastic  region  back  to  the  line 


<so 
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O 
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V.  s. 


s-/(j 


2k 


r  -  b  =  CgCt  -  2.A).  In  ihe  elagtic  region  the  third  diagonal  row 
is  started  at  (9>7-)  rather  than  at  (10,8).  Since  the  changes  which 
are  made  along  the  edge  of  the  plastic  region  will  affect  eisplace- 
ments  at  later  times,  it  is  impprtant  that  the  shearing  strain  at 
the  edge  of  the  plastic  region  be  evaluated  f3!*equently  as  one  moves 

^  t 

down  the  diagraja-  In  order  to  avoid  useless  calculation. 

In  our  numerical  work  we  have  used  the  following  values  for 


the  constants  T/diich  appear  t 
Cl  =  114,910  ft/sec. 


^-3 


7^  =  B.Sh  X  10-  , 


“  11,21(0  ft/sec. 


gA  =  6.9ii, 

Cg/ci  =  0.751(,  b  =  1  ft, 

G  =  Ib/ft  sec^  ~  5l8  ft/sec, 

23'q(G  -  F)//'=  1.61i  X  10^  ft® /sec? 


The  foregoing  values  of  the  constants  are  not  appropriate  for  armor 
plate j  however,  the  qualitative  nature  of  the  waves  produced  does 
not  depend  on  the  numerical  values  of  the  constants.  The  wave 
velocities  are  both  about  25,  percent  too  small,.  The  stresses  that 
we  calculate  will  therefore  be  scanewhat  too  large.  From  t  =  0  to 
t  =  T,.1  X  10“^  sec  the  time' interval  A  used  is  10“5  sec-  From 
t  =  1,1  X  10~^  to  t  =.3.3  X  10“^  sec  the  interval  used  is  2x  10“5  sec 
The  large  time  interval  is  chosen  so  that  the  change  in  the  dis~ 
placement  in  going  frcm  one  point  in  a  diagonal  row  to  the  next  is 
never  more  than  1 6  percent  of  the  displacement.  The  average  change 
is  much  less,  amounting  to  5  percent  of  the  displacement.  In 

'x 

Figs.  1(,  5, .and  6  we  h^ve  plotted  the  displacements  and  the  shearing 
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stresses  as  a  function  of  the  radius  at  times  t=  1.1  x10“^  seG_,  ■ 
t  =2.3  x10“^  sec  and  t  =  3.3><  10“^  sec.  figures  h,  5  and  6  show  that 
the  shearing  strain  decreases  as  the  radius  r  increases.  Initially 
there  is  a  sizeable  discontinuity  in  the  shearing  strain  at  the  edge 
of  the  plastic  region,  but  this  becomes  less  marked  as  time  passes. 

5.  Future  program 

From  what  has  been  said  it  is  clear  thtit  the  investigation  of 
the  propagation  of  plastic  deformation  in  cylinders  is  by  no  means 
complete.  In  this  section  we  shall  mention  a  fevi  points  that  should 
be  investigated. 

The  derivation  of  the  Yjave  equations  in  this  report  indicates 
that  there  should  be  a  minimum  wave  velocity  for  cylindrical  waves 
of  deformation.  This  miniraum  wave  velocity  depends  chiefly  on  the 
compressibility  and  the  density  of  the  material.  In  the  case  of 
the  propagation  of  plastic  deformation  in  a  vlre  no  such  minimum 
velocity  is  found.  It  would  be  interesting  to  determine  whether 
such  a  minimum  -wave  velocity  can  be  found  experimentally. 

In  the  penetration  of  armor  plate  the  strains  near  the  axis 
of  the  shell  are  not  infinitesmal.  Therefore,  the  effect  of  these 
finite  strains  on  the  propagation  of  plastic  deformation  near  the 
axis  of  the  expanding  cylinder  should  be  considered  theoretically. 

A  solution  of  the  problem  in  the  form  of  a  series  or  some 
approximate  function  w'ould  be  of  much  use  since  one  could  easily 
determine  the  result  of  a  change  in  the  physical  constants  of  the 
material  under  consideration.  For  this  reason  it  would  be  well  to 
attempt  an  analytical  solution  of  the  case  wrhere  the  cylinder  contains 
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both  elastic  and  plastic  'waves.  The  approximate  solution  in  the 
case  where  only  elastic  -mves  are  present  should  furnish  useful 
information  for  this  work. 

.  The  experience  gained  in  working  -vw-th  a  Wo-dimensional  model 
of  the  penetration  problem  might  conceivably  point  the  -way  to  a 
numerical  solution  of  the  three-dimensional  problem  of  the  pene¬ 
tration  of  armor  plate. 

In  the  n'umerical  solution  given  in  this  report  the  displace¬ 
ment  at  the  surface  of  the  expanding  cylinder  can  be  taken  to  be 
any  desired  function  of  the  time.  Thus,  one  could  investigate 
the  "secondary  ’.mve"^^  which  is  sent  out  when  the  displacement  at 
the  surface  of  the  expanding  cylinder  reaches  its  maximum  value. 

To  do  this  one  -would  simply  carry  through  a  numerical  solution 
using  the  wave  equations,  Eqs.  (U),  (8)  and  (13)^  together  vdth 
the  boundary  condition, 

u(b)  =  0,  for  t  <  0; 

.u(b)  =  f^i"  0  <  t  <  T; 

u(b}  =  VqT,  for  T  <  t. 

It  is  obvious,  of  course,  that  more  detailed  ntawerical  calculations 

should  be  carried  out  using  constants  appropriate  for  armor  plate. 
How'ever,  it  would  seem  that  such  numerical  calculations  should  be 

»• 

made  after  one  knows  how  to  treat  the  finite  strains  near  the  shell. 
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UNCLASSIFIED 


GLOSSARY 

b  The  initial  radius  of  the  cylindrical  hole  in  a  thick 

plate  of  infinite  extent  (Secs.  3  a-nd  [|) . 

Cl  ’Yave  velocity  in  the  elastic  region. 

Cg  Yilave  velocity  in  the  plastic  region. 

D  Slope  of  the  stress-strain  curve  in  the  plastic  region. 

E  Young’s  modulus. 

F  Modulus  in  plastic  region  corresponding  to  shear  modulus 

in  elastic  regionj  F  =  d/2(l  +  c* ) , 

a<s  :  o  , 

f  Slope  of  stress -^strain  curve,  soraa  function  of  the  largest 

shearing  strain. 

f '  Derivative  of  function  f  vd.th  respect  to  its  arguinent. 

G  Modulus  of  rigidityj  G  =  E/2(1  +  v  ) ',  also  referred  to  as 

the  shear  m-odulus . 

K  Modulus  of  volijrae  expansion. 

r,G  Plane  polar  coordinates. 

T  Time  during  vhich  expansion  of  hole  takes  place  (Sec.  3). 

t  Tine . 

u  Displacement. 

Vq  Radial  velocity  of  material  a.t  the  surface  of  the  expand¬ 

ing  cylindrical  hole  (Sec.  3)* 

Defined  by  2oc  =  1  -  d/3K 

A  Ti^ie  interval. 

Shearing  strains. 

) :)  Shearing  strain  virhen  material  begins  to  deform  plastically. 

)x  Shea,ring  strain  at  point  of  decreasing  strains. 

£r  ^0  ^2  Tensile  strains  in  the  directions  v,  9  a.nd 

tv.  ^00  jto;  Strains  present  vj'hen  the  material  begins  to  deform  plasti- 
”  o  0  cally. 
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j 


Strains  present  at  the  time  when  the  stresses  in  the 
vicinity  of  the  cylindrical  hole  begin  to  decrease. 


X 


V 


The  yield  stress  obtained  in  tensile  tests  made  on 
long  thin  wires. 

Poisson  ratio. 


Density  of  the  material. 


Tensile  stresses  in  the  directions  r  and  ©. 


S?  required  to  initiate  plastic  deformation. 

(Tj.  j<rQ  ,(1^  Stresses  present  at  the  time  when  the  stresses  in 
the  vicinity  of  the  cylindrical  hole  begin  to  de¬ 
crease. 
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In  the  present  paper  it  is  assumed,  in  accord  with  Bethe’s  model  for  ixmor  penetration, 
that  the  plastic  deformation  produced  during  armor  penetration  is  similar. tp  that  pro* 
duced  by  an  expanding  cylinder.  In  the  first  portion  of  the  paper,  the  relations  l^tween 
the  stresses  and  the  strains  are  considered,  and  the  wave  equations  which  govern  the 
motion  of  the  material  when  it  is  rapidly  deformed  are  derived.  It  is  found  that  in  the 
case  of  the  thick  plate,  an  elastic  wave  diverges  radially  from  the  expandii^  cylinder, 
that  this  is  followed  by  a  plastic  wave,  and  that  the  elastic  and  plastic  wave  velocities 
do  not  differ  very  much.  In  the  second  portion  of  the  paper,  an  approximate  e:qpression 
for  the  displacement  is  obtained  for  the  case  where  the  deformaUon  is  elastic.  In  the 
last  section  of  the  paper,  a  method  of  numerically  integratii^  the  wave  equations  is  given. 
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